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Abstract. A new technique is used to derive the extinction and age of popula- 
tions in the interacting galaxies NGC5194 & NGC5195 from high S/N multicolor 
O^ ' photometric data. A new evolutionary scenario of the interaction is proposed. 
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^ ! 1. Introduction 

(N 

^ . Detailed investigation of populations in Local Group (LG) galaxies 

CN \ using stellar photometry has recently become very popular (Aparicio 

*^ ■ 1998). This is a direct way to study the populations and evolution 

vQ . of galaxies. However, even for LG objects it is difficult to reach the 

^^ ' required photometry depth. To reconstruct the star formation history 

On . of more distant galaxies a different techniques should be applied. Histor- 

^^ \ ically, the most popular method is an evolutionary population synthesis 

i-Q ' (Tinsley 1980; Arimoto & Yoshii 1986), which treats the entire galaxy 



Ph' 



as a single zone unit. Such an approach is probably satisfactorily appli- 



O \ cable only for ellipticals. However, it seems too simplified in the case of 

-4— » ' spiral disks. To improve the tools available for the evolutionary analysis 

C^ , of galaxies, a new method for multicolor mapping is introduced and 

applied for the interacting galaxy pair NGC5194/95 (hereafter, M51 
. system). 

r> \ M51 is a well studied nearby system of two overlapping galaxies seen 

C^ ■ nearly face-on, and therefore suitable for testing new methods and 2 — D 

mapping techniques. Since the interpretation of the results of optical 
photometry always suffer from the influence of dust residing in galaxies, 
it is important to know the extinction properties of the galaxies under 
consideration. For this purpose, overlapping galaxies are particularly 
suitable, allowing precise derivation of an attenuation law and the 
amount of the extinction in the foreground galaxy (White & Keel 
1992). Moreover, results of an extensive dynamical modeling of M51 
are available (Barnes 1998), and hence can be used for a comparison of 
the photometrically derived parameters with the ones determined by 
other independent methods. 
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Figure 1. CMDs of the galaxies NGC5194 & NGC5195. Galaxies are simply sepa- 
rated by declination. Reddening arrow is given according to the M51 attenuation 
law. 



2. Observations and Reductions 

Observations were carried out at Kiso Observatory (Japan) using tlie 
1.05 m Sclimidt telescope on February 17 & 18 and March 7, 1994. 
The CCD camera Tek 1024 (read-out noise 14 e~, scale 0.75 " /pixel, 
and field of view 12.5' x 12.5') was used. The open cluster M67 was 
observed every night as a standard field. More than 150 exposures of 
M51 of 30 — 1800 seconds duration and image quality of 3" — 5" were 
taken in total. The standard BVRc'Ic a-i^d narrow-band Ha (Aq = 
658 nm; SX = 2.5 nm) filters, together with Z (Aq = 518 nm; 6X = 
21 nm) filter of Vilnius system. 

Standard reductions were made with IRAF. The 2 — D aperture 
photometry (scanning the entire field with a constant step of ~ 4" and 
aperture of ~ 8") was performed using an AFO package (Bridzius &: 
Vansevicius 1998). 

Assuming the distance to the M51 system to be equal to 8.4 Mpc 
(Feldmeier, Ciardullo, & Jacoby 1997), each circular ~ 8" aperture 
corresponds to the region of ~ 330 pc in galaxies (hereafter, zone). In 
the color-magnitude diagram (CMD) (Fig. 1) zones of both galaxies are 
shown separately. For further discussion we use only precisely measured 
(cj < 0.'^03) zones with V < 23.'"5. 

It is worth noting the different morphologies and fine structures of 
both CMDs (Fig. 1). Most of the "streaming features" are parallel to 
the reddening arrow, and can be attributed to the dust infiuence. 



3. Extinction in M51 

The problem of dustiness of galaxies at different evolutionary stages 
is widely discussed recently. In particular, the dilemma of whether the 
spiral disks are optically thick or thin is debated (Davies & Burstein 
1995). To clarify the problem of spiral disk dustiness, we determine 
the amount of extinction and the attenuation curve in NGC5194 us- 
ing NGC5195 as a tracer. The facts supporting the validity of such a 
method are: a) NGC5195 is behind NGC5194; b) NGC5195 is a source 
bright enough to have a high S/N ratio even in the most obscured parts 
of the galaxy; c) NGC5195 is found to be a regular SO galaxy in the 
near-infrared (Thronson, Rubin, &: Ksir 1991), therefore it is possible 
to compare symmetrically located zones. 
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Figure 2. V (panel (a): solid line - clear (without dust signatures) region, dashed 
line - dusty region) Ik B — V (panel (b): no shading - clear region, shaded histogram - 
dusty region) distributions of zones used to derive the attenuation law and extinction 
in the disk of NGC5194. 

Table I. Properties of the extinction in NGC5194. 



Parameter 


NGC5194 


Milky Way 


Rv 


2.56 ±0.13 


3.10 


Ev-r/Eb~v 


0.67 ±0.06 


0.56 


Ev-i/Eb-v 


1.33 ±0.08 


1.25 


Eb-i/Eb-v 


2.41 ±0.13 


2.25 


Er-i/Eb-v 


0.66 ±0.09 


0.69 


Ez-v/Eb-v 


0.30 ±0.05 


0.28 



We measured two regions (symmetric with respect to the NGC5195 
center and bar), each consisting of 50 zones (~ 20" x 40") at an average 
distance of ~ 40" from the center. Visual inspection of the B image 
helped to locate the region relatively free of dust signatures, and then 
the corresponding dusty region was chosen on the opposite side of the 
center and bar. The V and B — V histograms of the investigated zones 
are shown in Figure 2. The parameters of the derived attenuation curve 
are given in Table I (definition of the extinction: Ay = Rv • Eb-v-, is 
used). The derived median value of the extinction for the interarm 
region at ~ 11 kpc distance from the center of NGC5194 is Ay = 
0.7 ±0.3. 

Numerous determinations of the amount of extinction in NGC5194 
have recently been obtained. Smith et al. (1990) derived Ab = 0.9 
using NGC5195 as a tracing source. Rix &: Rieke (1993) put a limit 
on the global optical depth of NGC5194 disk of ry < 1 - 2. Nakai k 
Kuno (1995) determined Ay = 1-16 — 4.38, based on a study of 37 
HI I regions in the disk of NGC5194. Beckman et al. (1996) concluded 
that optical depth in the V band, averaged out to a radius of three 
scale lengths from the center of NGC5194, is 0.7. These results mainly 
ascribe properties of the arm regions, however they do not contradict 
our conclusion that the disk of NGC5194 is opaque in the interarm 
regions even at a large distance from the center. 

However, such a conclusion is in contradiction with the one drawn 
by White, Keel, & Conselice (1996) from a study of nine overlapping 
galaxy systems. They derived the extinction Ab = 0.3 — 2.3 in the arm 
and Ab = 0.08 — 0.47 in the interarm regions. NGC3314 is the only 
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system in their sample which possesses significant extinction in the arm 
{Ab = 1-11 — 1.64) and in the interarm {Ab = 0.77 — 1.75) regions, 
which supports our conclusion that some of the spiral disks are opaque 
even in the interarm regions. 

The color excess ratios derived for the outer part of the disk are 
very close to those of the Milky Way (MW) extinction law. Taking 
into account the same conclusion reported by Panagia et al. (1996) for 
the central part of NGC5194, we assume that color excess ratios are 
constant throughout the disk of this galaxy. 



4. Abundance and effective age of zones 

To avoid the age-abundance degeneracy in determining the effective 
age of zones, global distribution of the abundance should be derived. 
Abundance analysis of HII regions (Diaz et al. 1991; Zaritsky, Kenni- 
cutt, & Huchra 1994) indicates the mean abundance to be above solar: 
12 + log{0/H) ~ 9.3, with a gradient equal to ~ —0.06 dex/kpc. Al- 
ternatively, Hill et al. (1997) explained the observed UV color gradient 
in terms of radial extinction gradient, rather than of abundance. 

We define an abundance parameter (independent of the extinction) 
as follows: QzvBV = {Z — V) — {B — V)-Ez-v/Eb-v- Qzvbv is a mea- 
sure of Mg (518 nm) spectral feature strength, which is a good abun- 
dance indicator (Barbuy 1994; Casuso et al. 1996). Plotting Qzvbv 
vs. radial distances from the galaxy centers, we do not see the abun- 
dance gradient except for the very central region inside the radius of 
~ 1.5 kpc. However, a significant difference of abundances between 
NGC5194 and NGC5195 is found. Therefore, deriving the ages of zones 
we assume no abundance gradient across the disk of both galaxies. 
The abundances Z=0.05 for NGC5194 and Z=0.02 for NGC5195 are 
adopted, respectively. Such an assumption does not significantly affect 
differential analysis of the age distribution of populations in the disk 
outside the inner Lindblad resonance (radius ~ 1.25 kpc) of NGC5194 
(Elmegreen, Elmegreen, & Seiden 1989), and outside the comparable 
radius in NGC5195. 

Let us introduce a parameter Qsy/ = (B — V) — {V—I)-Eb~^v/Ev-i 
as a measure of the effective age of each zone, which is again extinction- 
independent. Actually, both Q formulae are only correct in the case of 
a star residing behind the dust cloud (screen). The validity of dusty 
screen approximation in the case of face-on galaxies was shown by 
Gonzalez & Graham (1996). Diagram V vs. Qbvi constructed from our 
data sample, in contrast to Figure 1, shows no prominent "streaming 
features" , which can be attributed to the influence of extinction. All this 
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Figure 3. a) Qbvi vs. log(t) for solar (solid line) & 2.5 times solar (dashed line) 
abundances; b) Qbvi of NGC5194 (solid line) & NGC5195 (dashed line). 

Figure 4- A map of the population ages (Qbvi) in M51. Light gray corresponds to 
the youngest and dark gray to the oldest zones. 



indicates that Qbvi is a reddening-free parameter in the case of face-on 
galaxies. Figure 3a shows cahbration of Qbvi vs. age for single (age) 
stellar populations (SSPs) calculated from the model spectra library 
GISSEL96 (Bruzual &: Chariot 1993). Figure 3b provides histograms 
of the parameter Qbvi for both galaxies of the M51 system, which 
are simply separated by declination. A significant difference in age 
between the two galaxies is obvious even if the interfering influence of 
the Northern arm of NGC5194 projected onto NGC5195 is neglected. 
The derived ages of the bulk of population in NGC5194 k NGC5195 
are 390 it 70 Myr and 710 it 120 Myr, respectively. The gray-coded 
population age distribution in galaxies is shown in Figure 4. 

Recent dynamic modeling suggests two starbursts in M51, the first 
less than ~ 100 Myr ago and the second one ~ 400 Myr ago (see 
Barnes (1998) for discussion). We derive the age for a bulk of zones in 
NGC5194 to be ~ 390 Myr, in good agreement with the dynamical 
prediction. Zones in NGC5195 are ~ 320 Myr older. If one assumes 
this age difference to be a characteristic time between two close-up 
passages, the predicted age ~ 70 Myr for the last passage again is 
in good agreement with the dynamical modeling. Therefore, our study 
indicates that the evolutionary scenario of M51 has had at least three 
periods of activity. The first starburst occurred ~ 700 Myr ago in 
NGC5195 during a close-up passage with NGC5194, when, probably, 
NGC5195 was transformed to an SO galaxy. The next starburst took 
place after ^^ 300 Myr, when the bulk of the population in NGC5194 
was formed, and finally the last one occurred again after ~ 300 Myr 
when the present grand-design appearance of NGC5194 was "touched 
up". 



5. Conclusions 

Main conclusions are the following: a) the attenuation law determined 
in NGC5194 is close to MW standard extinction law; b) the opaque 
interarm region at distance ~ 11 kpc from the center of NGC5194 
supports optically thick models of spiral disks; c) the derived effective 
age ~ 400 Myr for the bulk of NGC5194 zones agrees well with the 
predictions of recent dynamic modeling of M51; d) the effective age 
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~ 700 Alyr derived for the bulk of zones in NGC5195 implies a scenario 
with three periods of starburst activity. However, more sophisticated 
radiative transfer models and realistic star-dust distributions as well 
as 2 — D models of the evolution of spiral disks, instead of simple SSP 
approach, are crucial for the interpretation of the mapping data. 
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